Digital Coherent and Dissipative Quantum Simulations with Rydberg Atoms 
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Following Feynman and as elaborated on by Lloyd, a universal quantum simulator (QS) is a 
controlled quantum device which reproduces the dynamics of any other many particle quantum 
system with short range interactions. This dynamics can refer to both coherent Hamiltonian and 
dissipative open system evolution. We investigate how laser excited Rydberg atoms in large spacing 
optical or magnetic lattices can provide an efficient implementation of a universal QS for spin 
models involving (high order) n-body interactions. This includes the simulation of Hamiltonians 
of exotic spin models involving n-particle constraints such as the Kitaev toric code, color code, 
and lattice gauge theories with spin liquid phases. In addition, it provides the ingredients for 
dissipative preparation of entangled states based on engineering n-particle reservoir couplings. The 
key basic building blocks of our architecture are efficient and high-fidelity n-qubit entangling gates 
via auxiliary Rydberg atoms, including a possible dissipative time step via optical pumping. This 
allows to mimic the time evolution of the system by a sequence of fast, parallel and high-fidelity 
n-particle coherent and dissipative Rydberg gates. 
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I. INTRODUCTION 

Laser excited Rydberg atoms G10, [1,0, H, 0,0] stored 
in large spacing optical lattices |g] or magnetic trap ar- 
rays [9(| offer unique possibilities for implementing scal- 
able quantum information processors. In such a setup 
single atoms can be loaded and kept effectively frozen 
at each lattice site, with long-lived atomic ground states 
representing qubits or effective spin degrees of freedom. 
Lattice spacings of the order of a few /im allow single 
site addressing with laser light, and thus individual ma- 
nipulation and readout of atomic spins. Exciting atoms 
with lasers to high-lying Rydberg states and exploit- 
ing the strong and long-range dipole-dipole or Van der 
Waals interactions between Rydberg states provides fast 
and addressable 2-qubit entangling operations or effec- 
tive spin-spin interactions; recent theoretical proposals 
have extended Rydberg-based protocols towards a single 
step, high-fidelity entanglement of a mesoscopic number 
of atoms [lfj. 1 1 11] . Remarkably, the basic building blocks 
behind such a setup have been demonstrated recently in 
the laboratory by several groups [H, [H] ■ 

Motivated by and building on these new experimental 
possibilities, we discuss below a Rydberg QS for many 
body spin models. As a key ingredient of our setup (see 
Fig. 1) we introduce additional auxiliary qubit atoms in 
the lattice, which play a two-fold role: First, they control 
and mediate effective n-body spin interactions among a 
subset of n system spins residing in their neighborhood 
in the lattice. In our scheme this is achieved efficiently 
making use of single-site addressability and a parallelized 
multi-qubit gate, which is based on a combination of 
strong and long-range Rydberg interactions and electro- 
magnetically induced transparency (EIT), as suggested 
recently in Ref. Second, the auxiliary atoms can 

be optically pumped, thereby providing a dissipative el- 



ement, which in combination with Rydberg interactions 
results in effective collective dissipative dynamics of a set 
of spins located in the vicinity of the auxiliary particle, 
which itself eventually factors out from the system spin 
dynamics. The resulting coherent and dissipative dynam- 
ics on the lattice can be represented by, and thus simu- 
lates a master equation p — —(i/H) [H,p]+£p [14], where 
the Hamiltonian H = J2 a H a is the sum of n-body in- 
teraction terms, involving a quasi-local collection of spins 
in the lattice. The Liouvillian term C — J2p ^{ c p) with 
V{c)p = cpc^ — 2 C cp — P2 C c in Lindblad form governs 
the dissipative time evolution, where the many-particle 
quantum jump operators cp involve products of spin op- 
erators in a given neighborhood. 

The actual dynamics of our system is performed as a 
stroboscopic sequence of coherent and dissipative oper- 
ations involving the auxiliary Rydberg atoms over time 
steps r, with the master equation emerging as a coarse 
grained description of the time evolution. For purely co- 
herent dynamics governed by the Hamiltonian, this is the 
familiar "digital QS" [Hi fig where for each time step 
the evolution operator is implemented via a Trotter ex- 
pansion e r lHT l h rj Y\ a e~ lHaT / h and a certain error as- 
sociated with the non-commutativity of the quasi-local 
interactions H a . The concept of stroboscopic time evo- 
lution is readily adapted to the dissipative case by inter- 
spersing coherent propagation and dissipative time steps 
e Cr w Y[p e D ( c ^ T ', providing an overall simulation of the 
master equation by sweeping over the whole lattice with 
our coherent and dissipative operations. Many of these 
steps can in principle be done in a highly parallel way, 
rendering the time for a simulation step independent on 
the system size. 

In our scheme the characteristic energy scale of the 
many-body interaction terms is essentially the same for 
two-body, four- or higher-order interaction terms, and 
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Figure 1: Setup of the system: a) Two internal states \A)i and 
\B)i give rise to an effective spin degree of freedom. These 
states are coupled to a Rydberg state \R)i in two-photon res- 
onance, establishing an EIT condition. On the other hand, 
the control atom has two internal states |0} c and |1} C . The 
state |l) c can be coherently excited to a Rydberg state \r) c 
with Rabi frequency O r , and can be optically pumped into 
the state |0) c for initializing the control qubit. b) For the 
toric code, the system atoms are located on the links of a 2D 
square lattice, with the control qubits in the centre of each 
plaquette for the interaction A p and on the sites of the lattice 
for the interaction B a . Setup required for the implementation 
of the color code (c), and the U(l) lattice gauge theory (d). 

mainly limited by the fast time-scale to perform the par- 
allel mesoscopic Rydberg gate operations. We note that 
this is in contrast to the familiar analog simulation of 
Hubbard and spin dynamics of atoms in optical lattices 

nana. 

where collisional interactions between atoms 
provide naturally two-body interactions, while higher or- 
der, small effective interactions and constraints are typi- 
cally derived with perturbative arguments [l9| . 

II. COHERENT AND DISSIPATIVE MANY 
BODY SPIN DYNAMICS 

Before proceeding with the concrete physical imple- 
mentation of our Rydberg QS, we find it convenient to 
discuss special spin models and master equations of in- 
terest, starting with an explicit example: Kitaev's toric 
code. We will discuss the realization of a more complex 
setup of a three-dimensional U(l) lattice gauge theory 
giving rise to a spin liquid phase in the last section. 

Kitaev's toric code is a paradigmatic, exactly solvable 
model, out of a large class of spin models, which have 
recently attracted a lot of interest in the context of stud- 
ies on topological order and quantum computation. It 
considers a two-dimensional setup, where the spins are 
located on the edges of a square lattice [20(|. The Hamil- 

tonian H = —E ( ^ A p + J2 S Ba) is a sum of mutu- 



ally commuting stabilizer operators A p — Yiiep a f an< ^ 
B s = Yli^s °f > which describe four-body interactions be- 
tween spins located around plaquettes (Ap) and vertices 
(B s ) of the square lattice (see Fig. [TJd). The ground 
state (manifold) of the Hamiltonian consists of all states, 
which are simultaneous eigenstates of all stabilizer opera- 
tors A p and B s with eigenvalues +1, while the degeneracy 
of the ground state depends on the boundary conditions 
and topology of the setup. The excitations are charac- 
terized by a violation (eigenstates with eigenvalue — 1) 
of the low-energy constraint for each term A p ("magnetic 
charge") and B s ("electric charge") in the Hamiltonian, 
and exhibit mutual anyonic statistics. 

A dissipative dynamics which "cools" into the ground 
state manifold is provided by a Liouvillian with quantum 
jump operators, 

Cp=l*![l-A p ], c. = ~o*[l-B a ], (1) 

with i £ p and j £ s, which act on four spins lo- 
cated around plaquettes p and vertices s, respectively. 
The jump operators are readily understood as operators 
which "pump" from —1 into +1 eigenstates of the sta- 
bilizer operators: e.g. due to the projector (1 — A p )/2 
the jump operator c p will not act on states in the +1 
eigenspace of A p , while for —1 an arbitrary one of the 
four spins is flipped by erf. Similar arguments apply to 
c s . Efficient cooling is achieved by alternating the index 
i of the spin, which is flipped. The jump operators then 
give rise to a random walk of anyonic excitations on the 
lattice, and whenever two anyons of the same type meet 
they are annihilated, resulting in a cooling process, see 
Fig. [2j Note that by flipping, for example, the signs 
(1 — A p ) — > (1 + A v ) on selected plaquettes this dissipa- 
tive dynamics will prepare excited states. 

Our choice of the jump operator follows the idea of 
reservoir engineering of interacting many-body systems 
as discussed in Ref. [U, [22| . In contrast to alternative 
schemes for measurement based state preparation [23 |. 
here, the cooling is part of the time evolution of the sys- 
tem. These ideas can be readily generalized to more com- 
plex stabilizer states and to setups in higher dimensions, 
as in, e.g., the color code (see Fig. [Tfc) |2J |. 

III. IMPLEMENTATION OF A SINGLE TIME 
STEP 

We now turn to the physical implementation of the 
digital quantum simulation. The system and auxiliary 
atoms are stored in a deep optical lattice or magnetic trap 
arrays with one atom per lattice site, where the motion of 
the atoms is frozen and the remaining degree of freedom 
of the system and auxiliary atoms are effective spin— 1/2 
systems described by the two long-lived ground states 
\A)i and \B)i and |0) c and |l) c , respectively (see Fig. UK). 

We first discuss the elements of the digital QS for a 
small local setup, and present the extension to the macro- 
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Figure 2: Cooling of the toric code: a) A dissipative time 
step moves one anyonic excitation (red dot) on top of a sec- 
ond anyon sitting on a neighboring plaquette, annihilating 
each other and thus lowering the internal energy of the sys- 
tem. The anyon of different type (green dot) is unaffected, b) 
Numerical simulation of the cooling for JV lattice sites (peri- 
odic boundary conditions). Single trajectories for the anyon 
density n over time are shown as solid lines. Filled circles 
represent averages over 1000 trajectories. The initial state 
for the simulations is the fully polarized, experimentally eas- 
ily accessible state of all spins down. For perfect gates the 
energy of the system reaches the ground state energy in the 
long time limit, while for imperfect gates heating events can 
occur (blue solid line) and a finite density of anyons n remains 
present (blue circles). This finite anyon density corresponds 
to an effective temperature T off = — So/(fes log n). The pa- 
rameter quantifying the gate error was set to \Q\ = 0.1 (see 
App. [A] for details). 



scopic lattice system below. To be specific, we will focus 
on a single plaquette in the example of Kitaev's toric 
code outlined above. 

The implementation of the four-body spin interaction 
A p = Yli cf and the jump operator c p uses an auxiliary 
qubit located at the centre of the plaquette (see FigQJ)). 
The general approach then consists of three steps (see 
Fig. [3b) : (i) We first perform a gate sequence G which 
encodes the information whether the four spins are in a 
+1 or — 1 eigenstate of A p in the two internal states of the 
auxiliary atom, (ii) In a second step, we apply gate oper- 
ations, which depend on the internal state of the control 
qubit. Due to the previous mapping these manipulations 
of the control qubit are equivalent to manipulations on 
the subspaces with fixed eigenvalues of A p . (iii) Finally, 
the mapping G is reversed, and the control qubit is re- 
initialized incoherently in its internal state |0) by optical 
pumping. 

The mapping G is a sequence of three gate operations 



G = U c (ir/2)- 1 U g U c (ir/2), 



(2) 



where U c {tt/2) = exp(— wr<7 H /4) is a standard 7r/2-single 
qubit rotation of the control qubit and the parallelized 
many-body Rydberg gate [ll| takes the form (see Fig. UK 
for the required electronic level scheme and App. [El for 
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Figure 3: Single time step: a) The gate sequence G coherently 
maps the information, whether the system spins reside in any 
eigenstate [A, — ) (|A, +)) corresponding to the eigenvalue — 1 
(+1) of the many-body interaction A p onto the internal state 
|l)c (|0) c ) of the control qubit. b) After the mapping G, we 
apply gate operations, which depend on the internal state of 
the control qubit. Finally, the mapping G is reversed and 
the control qubit is incoherently reinitialized in state |0) c by 
optical pumping. At the end of the complete sequence the 
dynamics of the control qubit factors out. 



a brief summary) 

c/ 3 = |0)(0| c ®l + |l)(l| c 



(3) 



For the control qubit initially prepared in |0) c , the gate G 
coherently transfers the control qubit into the state |l) c 
(|0) c ) for any system state |A, — ) (|A, +)), with |A, ±) 
denoting the eigenstates of A p , i.e., A P \X, ±) = ±|A, ±), 
see Fig. d 

For the coherent time evolution, the application of a 
phase shift exp(i(f)a^) on the control qubit and the sub- 
sequent reversion of the gate, G _1 , implements the time 
evolution according to the many-body interaction A p , 
i.e., 



U iDt = exp (i(j)Ap) = G 1 exp (ificr*) G. 



(4) 



The control qubit returns to its initial state |0) c after the 
complete sequence and therefore effectively factors out 
from the dynamics of the system spins. For small phase 
imprints <j) <C 1, the mapping reduces to the standard 
equation for coherent time evolution 



dtp 



-E [-A p ,p]+o(4> 2 ). 



(5) 



The energy scale for the four-body interaction A p be- 
comes Eo = h(f>/T with r the time required for the im- 
plementation of a single time step. 

On the other hand, for the dissipative dynamics, we 
are interested in implementing the jump operator c p (see 
Eq. [1]). To this purpose, after the mapping G, we apply 
a controlled spin flip onto one of the four system spins, 



U Z M = |0)(0| C S)1 + |1)(1| C ®E 



(6) 
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with £ = exp(i#of). As desired, the sequence 
G~ 1 U Zii (0)G leaves the low energy sector |A, +) invari- 
ant, while - with a certain probability - it performs a 
controlled spin flip on the high energy sector |A, — ). How- 
ever, once a spin is flipped, the auxiliary qubit remains 
in the state |1} C , and optical pumping from |1} C to |0) c 
is required to re-initialize the system, guaranteeing that 
the control qubit again factors out from the system dy- 
namics. The optical pumping constitutes the dissipative 
element in the system and allows one to remove entropy 
in order to cool the system. The two qubit gate U ZA (0) 
is implemented in close analogy to the many-body Ryd- 
berg gate U g previously discussed. In one dissipative time 
step the density matrix p of the spin system is changed 
according to 



ApA ] + BpB ] 



with 



(7) 



(8) 



A = - [(A p + If + (A p - 1)E(A P + 1)] 

B = l[(A p -l)(A p + l) + (A p + l)E(A p -l)] 



For small phases 0, the change of the density matrix can 
be expressed by its time derivative and the mapping J7]) 
reduces to the Lindblad form 



d t p = k 



1 



vP c l ~ o { c I c pP + P c I c p} 



o(6*) (9) 



with the jump operators c p given in Eq. (TJ and the cool- 
ing rate k — 8 2 /t. Here, r denotes again the time scale 
for the implementation of a single time step. 

The above scheme for the implementation of the many- 
body interaction A p and the dissipative cooling with c p 
can be naturally extended to arbitrary many-body inter- 
actions between the system spins surrounding the con- 
trol atom, as e.g., the B p interaction terms in the above 
toric code. Gate operations on single system spins allow 
to transform of in of and of, while selecting only cer- 
tain spins to participate in the many-body gate via local 
addressability gives rise to the identity operator for the 
non-participating spins. Consequently, we immediately 
obtain the implementation of the general many-body in- 
teraction and jump operators 



-4, 



C0 



:Qi 



(10) 



with Wi,Qi E {1, erf , of , of }. Here, a and (3 stand for 
a collection of indices characterizing the position of the 
local interaction and the interaction type. Note that A a 
also includes single particle terms, as well as two-body 
interactions, while the dissipative term with Qi — 1 de- 
scribes simple dephasing instead of cooling. 



IV. TOOLBOX FOR DIGITAL QUANTUM 
SIMULATION 



Extending the analysis to a large lattice system with 
different, possibly non-commuting interaction terms in 
the Hamiltonian, i.e., H = ^2 a E a A a and dissipative 
dynamics described by a set of jump operators cp with 
damping rates Kp, provides a complete toolbox for the 
quantum simulation of many-body systems. Each term 
is characterized by a phase 4> a (Op) written during a single 
time step determining its coupling energy E a — Hcj> a /T 
and damping rate up — Op/r. For small phases cf) a <C 
1 and 6g <C 1, the sequential application of the gate 
operations for all interaction and damping terms reduces 
to the master equation of Lindblad form, 



d tP = -~[H,p] 



cppc' p 



\( 



epepp + pc\cp 



(11) 

The choice of the different phases during each time step 
allows for the control of the relative interaction strength 
of the different terms, as well as the simulation of inho- 
mogeneous and time dependent systems. 

The characteristic energy scale for the interactions E a 
and damping rates up are determined by the ratio be- 
tween the time scale r required to perform a single time 
step, and the phase difference <p a and 9p written during 
these time steps. It is important to stress that within our 
setup, the interactions are quasi- local and only influence 
the spins surrounding the control qubit. Consequently, 
the lattice system can be divided into a set of sublat- 
tices on which all gate operations that are needed for a 
single time step r, can be carried out in parallel. Then, 
the time scale for a single step r becomes independent 
on the system size and is determined by the product of 
the number z of such sublattices and the duration r s of 
all gate operations on one sublattice. In our setup, t s is 
mainly limited by the duration of the many-body Ryd- 
berg gate U g , which is on the order of a few hundred ns 
(see App. [B]for details). For the toric code discussed 
above, we have to apply the many-body gate twice for 
every interaction term (see Fig. [3]), and with z = 4, we 
obtain r ~ a few ps, resulting in characteristic energy 
scales and cooling rates of the order of several hundred 
kHz. It is a crucial aspect of this quantum simulation 
with Rydberg atoms that it can be performed fast and is 
compatible with current experimental time scales of cold 
atomic gases [25| . 

Finally, we would like to point out that imperfect 
gate operations provide in leading order small perturba- 
tions for the Hamiltonian dynamics and weak dissipative 
terms; see App. [A] for a detailed discussion and Fig. [2] for 
a numerical analysis of the induced errors. However, the 
thermodynamic properties and dynamical behaviour of a 
strongly interacting many-body system are in general ro- 
bust to small perturbations in the Hamiltonian and weak 
coupling to noise fluctuations. Consequently, small im- 
perfections in the implementation of the gate operations 
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are tolerable and still allow for the efficient quantum sim- 
ulation. 

An important aspect for the characterization of the 
final state is the measurement of correlation functions 
X = {A ai ...A an ), where A aj denote local, mutually 
commuting many-body observables. In our scheme, the 
observables A aj can be measured via the mapping G of 
the system information onto auxiliary qubits and their 
subsequent state selective detection. In analogy to noise 
correlation measurements in cold atomic gases [IE H3] 
the repeated measurement via such a detection scheme 
provides the full distribution function for the observables, 
and therefore allows to determine the correlation function 
X in the system. Consequently, in the above discussion of 
Kitaev's toric code, the necessary string operators char- 
acterizing topological order can be detected. 

V. LATTICE GAUGE THEORY 

In the first example, we discussed the implementation 
of the quantum simulator for the toric code, where the 
ground state is given by a stabilizer state, and the ex- 
tension to more complex stabilizer states is straightfor- 
ward. In the following, we will show that our approach 
can also be extended to systems with non-commuting 
terms in the Hamiltonian. As an example, we focus on 
a three-dimensional [/(l)-lattice gauge theory j28|, and 
show that dissipative ground state cooling can also be 
achieved in such complex models. Such models have 
attracted a lot of recent interest in the search for 'ex- 
otic' phases and spin liquids [IE H3, El, E2] ■ The three- 
dimensional setup consists of spins located on the links 
of a cubic lattice (see Fig. [TJi). The lattice structure 
for the spins can be viewed as a corner sharing lattice of 
octahedra with one site of the cubic lattice in the center 
of each octahedra. The Hamiltonian for the U{1) lattice 
gauge theory takes the form 

H = uY / (Sof- JJ2 B p + VN ™> (12) 

o p 

where the first term in the Hamiltonian defines a low 
energy sector consisting of allowed spin configurations 
with an equal number of up and down spins on each 
octahedron, i.e., spin configurations with vanishing to- 
tal spin S„ — J2i£o a i on eac h octahedron. The second 
term denotes a ring exchange interaction on each plaque- 
tte with B p = S+S^S+Si + S^S+S^S^; here Sf = 
[of ± ia\ ] /2 and the numbering is clockwise around the 
plaquette. This term flips a state with alternating up 
and down spins on a plaquette, i.e., | |, j, T, l) P — ► | J., T 
, J., T) P - The last term denotes the the so-called Rokhsar- 
Kivelson term, which counts the total number of flipable 
plaquettes 7V RK = J2 P Bp\ for a non-flipable plaquette 
vanishes. While the ring exchange interaction commutes 
with the spin constraint, ring exchange terms on neigh- 
boring plaquettes are non-commuting. At the Rokhsar- 
Kivelson point with J = V, the system becomes exactly 



solvable and it has been proposed that in the regime 
< V < J the ground state is determined by a spin liq- 
uid smoothly connected to the Rokhsar-Kivelson point 
[3ll |: the properties of this spin liquid are given by an 
artificial 'photon' mode, gapped excitations carrying an 
'electric' charge (violation of the constraint on an octa- 
hedron), which interact with a 1/r Coulomb potential 
mediated by the artificial photons, and gapped magnetic 
monopoles. 




Figure 4: Lattice gauge theory: a) Illustration of a dimer 
covering with three dimers (red links) meeting at each site 
of the cubic lattice. The front plaquette represents a fli- 
pable plaquette, which is transformed under the action of the 
ring-exchange B v into a different dimer covering, b) Numer- 
ical simulation for the cooling into the ground state at the 
Rokhsar-Kivelson point with E = for a system with 4 unit 
cells (12 spins). The cooling into the constraint on the octahe- 
dra follows in analogy to the cooling of the toric code via the 
diffusion and annihilation of 'electric charges' on the octahe- 
dra. The inset shows the cooling into the equal superposition 
of all dimer coverings starting from an initial state satisfying 
the constraint on all octahedra. c) Coherent time evolution 
from the Rokhsar-Kivelson point with a linear ramp of the 
Rokhsar-Kivelson term V(t) = J(l — tj/10h): the solid line 
denotes the exact ground state energy, while the dots repre- 
sent the digital time evolution during an adiabatic ramp for 
different phases <j) written during each time step. The differ- 
ence accounts for errors induced by the Trotter expansion due 
to the non-commutative terms in the Hamiltonian. 

In the following, we present the implementation of this 
Hamiltonian within our scheme for the digital quantum 
simulation and demonstrate that dissipative ground state 
cooling can be achieved at the Rokhsar-Kivelson point, 
from where the entire phase diagram is accessible. The 
control qubits reside in the center of each octahedron 
(on the lattice sites of the 3D cubic lattice) controlling 
the interaction on each octahedron, and in the center of 
each plaquette for the ring exchange interaction B p , see 
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Fig. [TJ Then, the coherent time evolution of the Hamil- 
tonian (fT2| can be implemented in analogy to the above 
discussion by noting that the ring exchange interaction 
Bp and N RK can be written as a sum of four-body in- 
teractions of the form lfT0|) . while the constraint on the 
octahedra is an Ising interaction, see App. Next, 
we discuss the jump operators for the dissipative ground 
state preparation. The cooling into the subspace with an 
equal number of up and down spins on each octahedron 
is obtained by the jump operator 



(13) 



where the product is carried out over the six spins located 
on the corners of the octahedron (see Fig. [TJi). The "in- 
terrogation" part 1 — TJ ■ ex P(*f a j) of the jump operator 
vanishes if applied to any state with three up and three 
down spins, while in all other cases a spin is flipped. Then 
the cooling follows in analogy to the cooling in the toric 
code by the diffusion of the 'electric' charges. Identifying 
each spin up with a 'dimer' on the link, all states satis- 
fying the constraints on the octahedra can be viewed as 
a dimer covering with three dimers meeting at each site 
of the cubic lattice, see Fig. Within this description, 
the ground state at the Rokhsar-Kivelson point is given 
by the condensation of the dimer coverings [12], i.e., the 
equal weight superposition of all dimer coverings. The 
condensation of the dimer coverings is then achieved by 
the jump operator 



(14) 



The action of this jump operator can be understood by 
starting with a dimer covering: For a non-flipable plaque- 
tte p (0 eigenstate of B p ), the dimer covering is a dark 
state of the jump operator c p , while for a flipable plaque- 
tte p (±1 eigenstate of B p ) the state is cooled into the 
equal weight superposition of the original dimer covering 
and the dimer covering obtained by flipping the plaque- 
tte (i.e., the +1 eigenstate). Consequently, the system 
is cooled into the dark state which is the equal super- 
position of all dimer coverings, which can be reached by 
flipping different plaquettes. The cooling of these jump 
operators is demonstrated via a numerical simulation for 
a small system of 4 unit cells, see Fig. HJd. 

The implementation of the digital quantum simula- 
tions provides full control on the spatial and temporal in- 
teraction strengths. Therefore, there are two possibilities 
to analyze the phase diagram for arbitrary interaction 
strengths: (i) The possibility to vary the different cou- 
pling strengths in time allows to adiabatically vary the 
different interaction terms from the dissipatively cooled 
ground state at the Rokhsar-Kivelson point; the adia- 
batic preparation using the Trotter expansion is shown 
in Fig. 3h. (ii) On the other hand, the spatial control 
of the coupling parameters allows us to divide the lat- 
tice into a system and a bath. The ground state of the 



bath is given by the Rokhsar-Kivelson state, which can 
be continuously cooled via the dissipative terms, while 
the system part is sympathetically cooled due to its con- 
tact with the bath; in analogy to the cooling well known 
in condensed matter systems. 
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Appendix A: GATE ERRORS 

In the following, we discuss the influence of a gate error 
onto the dynamics of the system. For simplicity, we illus- 
trate the general behaviour for an error in the many-body 
gate U g for the coherent time evolution of the many-body 
interaction A p . The imperfect many-body gate operation 
can be written 



U g = \0)(0\ c ®e i * c * + \1)(1\ C 



■A r 



(Al) 



where the perfect gate U g is recovered for Q — > and the 
operator Q = acts on the system spins surrounding 
the control atom. This form of the error is motivated 
by the specific implementation of the gate [ll|; however, 
it can be seen that different errors in the many-body 
and single particle gates will lead to similar phenomena. 
Note that we have introduced <f> into the definition of the 
error in order to perform a consistent expansion. For the 
coherent time evolution, the imperfect gate gives rise to 
a finite amplitude for the control qubit to end up in the 
state |l) c . Consequently, optical pumping of the control 
qubit is required to reinitialize the system. Then the gate 
operations on a single plaquette give rise to the mapping 
of the density matrix onto 



P - 



CpC^ + DpD^ 



(A2) 



with (9 = e w ) 



C = - [(e 2 + l)cos</>+(eA p + A p e)isin0](A3) 

« exp[ic/>(A p + Q)}-^Q 2 (A4) 

D = i [(e 2 -l)cos</>+(eAp-A p e)isin0](A5) 
« -i<f>Q. (A6) 

The last equations hold with an accuracy up to third or- 
der in the small parameter 4>. Consequently, the optical 
pumping has no influence in leading order, and the sys- 
tem is well described by a Hamiltonian evolution with 
the modified Hamiltonian H = —{h/T)<j) [Ap + Q]. The 



7 



characteristic energy scale of the correction is again given 
by htfi/r, and consequently describes a small perturbation 
if \Q\ <C 1. In second order expansion in the small pa- 
rameter 4>, the mapping of the density matrix reduces 
to 



p-> p-i<p[h, p] 



[h, [h,p]] 



(2QpQ-{Q 2 ,p}) 



2 L "' L "' rjJ ' 2 

(A7) 

with h = — [A p + Q] . The first terms on the right hand 
side describe the coherent evolution of the system with 
the evolution operator exp(—i(j>h) consistently expanded 
up to second order, while the last term takes the standard 
Lindblad form for a dissipative coupling with the jump 
operator c e = Q describing a dephasing with the rate 
n e = 4> 2 jr. The dissipative terms appear only in second 
order and are therefore strongly suppressed for nearly 
perfect gates with \Q\ <C 1. 



the required lattice spacing of the optical lattice, which 
is essentially given by the distance where the dimension- 
less interaction V — AAC^/^^ 6 ) equals one, is given by 
r = 1.4 pm. This lattice spacing can be achieved using 
standard techniques such as tuning the angle between the 
lattice laser beams, and it is also larger than the LeRoy 
radius, i.e. the interaction between the Rydberg atoms 
is well described by a van der Waals interaction. 



Appendix C: QUANTUM GATES FOR THE (7(1) 
LATTICE GAUGE THEORY 



The Hamiltonian giving rise to the constraint for the 
spins on the octahedra can be expressed as a sum of Ising 
interactions, 



Appendix B: MESOSCOPIC RYDBERG GATE 

In the following we briefly summarise the main prop- 
erties and requirements of the many-body Rydberg gate 
U g introduced in Ref . . The internal level structure of 
the control atom and the surrounding ensemble atoms is 
depicted in Fig. [TJ The underlying physical mechanism 
of the gate operation ([3]) is a conditional Raman transfer 
of all ensemble atoms between their logical internal states 
|^4) and \B), which - depending on the internal state |0) 
or |1) of the control qubit - is either inhibited or enabled. 
The gate is realised by the following three laser pulses: (i) 
A first state selective 7r-pulse acting on the control atom 
changes the ground state |1) into the Rydberg state \r). 
(ii) During the whole gate operation, a strong coupling 
laser of Rabi frequency fl c constantly acts on all ensem- 
ble atoms and off-resonantly couples the Rydberg level 
\R) to the intermediate level \P) with a detuning A. Its 
frequency is chosen such that is in two-photon resonance 
with the two Raman laser beams of Rabi frequency ft p 
(see Fig. [J) , thereby establishing a condition known as 
electromagnetically induced transparency (EIT) (34[. In 
consequence, when the Raman laser pulses are applied 
- and provided the control atom resides in state |0) - 
this EIT condition effectively inhibits the coupling of the 
Raman lasers to the intermediate state \P) and the Ra- 
man transfer is blocked. In case the control atom was 
excited to the Rydberg state |r) in step (i), the large 
Rydberg-Rydberg interaction energy shift (dipole block- 
ade) lifts the EIT condition for the ensemble atoms and 
thus the Raman transfer takes place, (iii) Finally, the 
control atom is transfered from state \r) back to |1) via a 
second 7r-pulse. The total time T gate required for the gate 
is mainly limited by the duration of the Raman pulse, 
resulting in T gate ~ 16ttA/(302). For 87 Rb experimen- 
tally realistic parameters are a principal quantum num- 



oj + const, 



(CI) 



which allow for an efficient implementation using the gen- 
eral toolbox for quantum simulation. The implementa- 
tion for the jump operators for the constraint is obtained 
in analogy to the general jump operators with the many- 



body gate U g replaced by the gate 



|l)(l|c® 



fli exp(z^CT,f ). On the other hand, the ring exchange 
interaction can be written as a sum of commuting four- 
body interactions 



B,, 



8 ^ p 



(C2) 



Likewise, the Rokhsar-Kivelson term can be decomposed 
into 



3=1 
i J z z 

+a 1 a 2 a 3 a i 
+a 1 a 2 a 3 a i 



„zzJ^Ji 

a 1 a 2 a 3 a 4: 



„Q„0„z„z 

a 1 a 2 a 3 a 4 



a x a 2 a 3 a 4 
°i°*2<>i°i), (C3) 



where of is the identity matrix. Consequently, the co- 
herent time evolution follows again from the general tool- 
box, while the jump operators for cooling into the ground 
state at the Rokhsar-Kivelson point effectively cool into 
the zero eigenvalue eigenstate of the operators 



tally realistic parameters are a principal quantum num- ^ 1 16 1 

ber of n = 43 for the \R) state, % = 2tt x 100 MHz, -[1 - B P ]B P = — Cp } = ^ 



A = 2-7T x 1.2 GHz, and fi c = 2n x 1 GHz, resulting in 
a gate speed of T ga te ~ 320 ns. For these parameters 



3=1 i'=l 

(C4) 
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This can be achieved by replacing the gate U g with 
U B = |0)(0| c (8)l + |l)(l| c (g)exp[i|(l-B p )S p 

16 

= Y[ U c (tt/2)- 1 U j U c (tt/2) exp ^/32ct c 2 ) 

3=1 

U c {K/2)- 1 U 3 U c (n/2), (C5) 



with C/j = |0){0| C ®1 + |1)(1| C ®C^ ) . This gate operation 
leaves states with eigenvalue 0, +1 of B p invariant, while 
the —1 eigenvalue picks up a phase of 7T. It can be im- 
plemented as a product of many-body gates which derive 
directly from the standard gate U g with the combination 
of spin rotations. 
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